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We consider the collisions of plane gravitational and electromagnetic waves with 
distinct wavefronts and of arbitrary polarizations in a Minkowski background. We 
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(3JT), first present a new, completely geometric formulation of the characteristic initial 

value problem for solutions in the wave interaction region for which initial data are 
those associated with the approaching waves. We present also a general approach to 
the solution of this problem which enables us in principle to construct solutions in 
terms of the specified initial data. This is achieved by re-formulating the nonlinear 
dynamical equations for waves in terms of an associated linear problem on the 
spectral plane. A system of linear integral "evolution" equations which solve this 
spectral problem for specified initial data is constructed. It is then demonstrated 
explicitly how various colliding plane wave space-times can be constructed from 
given characteristic initial data. 
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1 Introduction 



The collisions of plane gravitational or gravitational and electromagnetic waves in general 
relativity have been studied extensively over recent decades by many authors. These stud- 
ies have revealed many interesting properties concerning the nonlinear interaction between 
such waves. These relate to both the physical and geometrical features of space-times with 
colliding plane waves. The great majority of these properties have been deduced from de- 
tailed studies of particular solutions of Einstein's equations, which were derived by some 
very special particular methods and tricks. However, a general mathematical solution 
of this problem, in which solutions in the wave interaction region are constructed from 
"input data" (the parameters characterizing the waves before their collision) has not been 
presented until now (although reported briefly in pQ). 

In this paper, we present a general scheme for a solution of this problem for gravita- 
tional or gravitational and electromagnetic waves of arbitrary polarizations propagating 
with distinct wavefronts and colliding in a Minkowski background. This method is based 
on the well known integrability properties of the hyperbolic reductions of Einstein's vac- 
uum and electrovacuum field equations with a two-dimensional Abelian space-time isome- 
try group. Mathematical tools for the solution of the characteristic initial value problems 
for these equations have been developed in 0. However, this approach has had to be 
adapted (with some necessary generalizations) for colliding plane waves in a Minkowski 
background. 

In the remaining part of this introduction, we recall some elements of the theory 
of integrable hyperbolic reductions of the Einstein and Einstein-Maxwell equations and 
give some references to preceding results. In the subsequent sections, we present a detail 
description of our new construction, the corresponding solution generation procedure, and 
we demonstrate how the procedure can be implemented. This demonstration includes the 
derivation of a previously unknown family of solutions. 

Reviewing what is known about colliding plane waves, one could initially recall an 
appreciable number of known solutions of Einstein's equations which can be interpreted 
as the outcomes of a collision of plane waves. A large number of such solutions and 
their physical and geometrical interpretations can be found in This list may be 
extended by further publications appearing since 1991. In particular, infinite hierarchies 
of exact vacuum and electrovacuum solutions with an arbitrary number of free parameters 
have recently been found j3], and many of these are of the type that is appropriate for 
colliding plane waves. However, all of these colliding plane wave solutions have been 
found using an "inverse" approach in which a formal solution in the interaction region is 
directly obtained and the initial parameters of the corresponding approaching waves were 
calculated afterwards. 

These solutions describe pairs of plane pure gravitational or mixed gravitational and 
electromagnetic waves with distinct wavefronts which approach each other from opposite 
spatial directions on the Minkowski background. The problem of constructing the corre- 
sponding solution in the wave interaction region is a well formulated characteristic initial 
value problem. (See |3| and the sections below for a detailed formulation of this problem 
including the general properties of approaching waves, typical matching conditions and 
the structure of the governing symmetry reduced field equations.) 

In the important particular case in which the colliding waves are purely gravitational 
and possess constant and aligned polarizations, Einstein's equations in the interaction 
region can be reduced to a linear Euler-Poisson-Darboux equation. In this case, the 
solution can be constructed from the characteristic initial data using the generalized 
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version of Abel's transform [S], However, if the polarizations of the approaching 
gravitational waves are not constant and aligned, or in the presence of electromagnetic 
waves, the governing Einstein or Einstein-Maxwell field equations are essentially nonlinear 
and this makes the problem much more complicated. 

Fortunately, these equations (which can be represented conveniently in the form of 
hyperbolic Ernst equations) have been found to be completely integrable. In fact, the 
integrability of the hyperbolic vacuum and electrovacuum Ernst equations facilitates a 
number of solution generating techniques such as the generation of vacuum [Zj and elec- 
trovacuum [H] solitons, Backlund or the symmetry transformations However, 
no techniques have been developed which provide algorithms for explicitly generating so- 
lutions from initial data. Nonetheless, this integrability gives rise also to some powerful 
methods for a general analysis of the governing field equations. In particular, Hauser and 
Ernst analyzed the characteristic initial value problem for the vacuum hyperbolic Ernst 
equation [T2]. They generalized their group-theoretical approach (which had been devel- 
oped earlier for stationary axisymmetric fields) and constructed a homogeneous Hilbert 
problem with corresponding matrix linear integral equations. Many global aspects of the 
characteristic initial value for vacuum fields, including the existence and uniqueness of 
solutions and a detailed proof of the Geroch conjecture, have been elaborated in this way 
However, this approach does not immediately lead to any effective methods for the 
solution of the corresponding characteristic initial value problem. 

Other schemes for the solution of characteristic initial value problems for the vacuum 
and electrovavuum Ernst equations ^1] or for the a-model form of the vacuum equations 
[To"] have also been considered. However, in these cases, characteristic initial value prob- 
lems were analyzed only for initial data which are regular functions of the geometrically 
defined coordinates (HJ) below. As will be discussed below, this restriction does not cover 
the collisions of plane waves with distinct wavefronts in a Minkowski background. In such 
cases, physically relevant initial data, being regular on the wavefronts in the local frame, 
must possess some singularities as functions of the coordinates adopted. 

The solution of the characteristic initial value problem for colliding plane waves that 
will be presented here originates from another approach to the analysis of the structure 
of integrable reductions of the Einstein and Einstein-Maxwell field equations that was 
developed in [TO] — - In this approach, every solution can be characterized by a set 
of functions of an auxiliary (spectral) parameter. These functions are interpreted as the 
monodromy data on the spectral plane of the fundamental solution of an associated linear 
system with a spectral parameter. These monodromy data are nonevolving (i.e. coordi- 
nate independent) and, generally, can be chosen arbitrarily or specified in accordance with 
the properties of the solution being sought. In particular, these data can be determined 
(at least in principle) from the initial or boundary data. In this scheme, the solution of 
the initial or boundary value problem is determined by the solution of some linear singular 
integral equations whose scalar kernel is constructed using these (specified) monodromy 
data. In this original form of the "monodromy transform" approach, a conjecture of local 
analyticity of solutions near the point of normalization was used very essentially. However, 
for waves with distinct wavefronts, the solutions of the linear equations are normalized 
at the point at which the waves collide and, for this situation, these are nonanalytic at 
this point. This violation of analyticity at the point of normalization gives rise to some 
difficulties with the construction of the master singular integral equation. 

Our resolution of this difficulty is achieved by the introduction of some "dressing" 
or "scattering" matrices (see [2]) into the structure of the fundamental solution of the 
associated linear system. The problem can then be overcome because, even in the non- 
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analytical case, these matrices possess much more convenient analytical properties. For 
locally analytic solutions, this construction is completely equivalent to the previous ap- 
proach. However, the new construction introduces some interesting features. In partic- 
ular, the monodromy data which characterize the scattering matrices are not conserved 
quantities. They evolve, and are therefore referred to as "dynamical monodromy data". 
Moreover, the evolution of these data is determined completely by the behaviour of fields 
on the characteristics which cross at the point of normalization. Again, the usual nonlin- 
ear partial differential form of the field equations are replaced by a pair of linear integral 
equations which determine the solution uniquely. With this new construction, the scalar 
kernels of these new quasi-Fredholm equations are built from this new kind of monodromy 
data. These linear integral equations then essentially possess a different, simpler struc- 
ture. Moreover, they are much better adapted to the context of a characteristic initial 
value problem. 

2 Space-time geometry for colliding plane waves 

Symmetry conjecture. In the situation under consideration, each of the approaching 
plane waves possesses two commuting space-like Killing vector fields with noncompact 
orbits, and we expect this symmetry to be global and retained during the collision and 
subsequent interaction of the waves. Thus, throughout the space-time, the metric and 
all nonmetric field components can be considered locally to be functions of two null 
coordinates, say u and v, while the two other space-like coordinates, say y and z (which 
we assume to be defined globally) are completely ignorable. An additional simplification 
of the metric and electromagnetic field components arises if we extend the conjecture of 
y- and z- independence from the metric and Maxwell tensor components to include their 
related potentials. 

Conjectured structure of Einstein Maxwell fields. For both colliding plane grav- 
itational waves in vacuum and a coupled system of colliding plane gravitational and elec- 
tromagnetic waves, the components of the metric and the 1-form of the electromagnetic 
vector potential can be considered globally in the forms 

ds 2 = g llu dx fJ, dx u + g ab dx a dx b , A = A b dx b (1) 

where //, is, . . . = 0,1 with a; , a; 1 = u, v, and a,b,... = 2,3 with x 2 ,x 3 = y,z. The 
orbit space metric g^ v , the metric on the orbits g ab and the components A b of the vector 
electromagnetic potential are functions of x^ only. 

General and geometrically defined coordinates. As is well known, the local coor- 
dinates x^ can be specified to reduce g^ v to a conformally flat form: 

g^ v dx^dx u = f(u,v)dudv. (2) 

It is clear that these null coordinates are defined up to an arbitrary transformation of the 
form 

u — u{u) and v = v(v). (3) 

For the analysis of the reduced Einstein-Maxwell equations and their integrability 
properties in the wave interaction region, it is convenient to use also some geometrically 



defined null coordinates £ and r\ which can be introduced as follows: 



i = (3 + a 
r) = (3 — a 



where 



a 



y/det \\g, 



ab 



and 



d v (3 = —d v a. 



(4) 



The function a(u,v) is the area measure on the orbits of the isometry group and the 
function j3{u,v) is then determined (up to an additive real constant). (It may be recalled 
that the Einstein-Maxwell equations imply that the function a(u,v) is a "harmonic" 
function in a sense that it satisfies the equation d u d v a = which also provides the 
integrability condition for the above /3-equations to be satisfied.) 

Parametrization of the metric components. It is convenient to parametrize the 
metric on the orbits by three scalar functions a(u, v) > 0, H(u,v) > and Q(u, v) as 



H 



HQ 



9ab ~ l HQ + - 

H, 



(5) 



Newman— Penrose tetrad and scalars. With the metric (0) it is convenient to adopt 
the Newman-Penrose null tetrad whose vector components with respect to coordinates 
(it, v, y, z) are: 



, 0,0,0} 



{0,1,0,0} 



n' 



|{0, 1,0,0} 



{1,0,0,0} 



m. 



|{o,o,i,n + ^} 



i H , „ ia 



where j = 0, 1, 2, 3. The expressions for the standard projections of the self-dual parts of 
the Weyl tensor (\l/ , \1> 2 , ^/ 3 , and of the Maxwell tensor (0 O , 1; 2 ) then take 
the forms 



(d v + ia~ 1 Hd v n) 



a 



#i = 



jd v (Q + -) 



VfH 



H 2 



2fa* 



d v (n + —)d u (n-f)- 



d u ad v a 



*3 = 



VfH 



(6) 



-(d u -ia- l Hd u tt) 



a 



jd u (Q - — ) 



where $ is the Ernst scalar electromagnetic potential to be defined below. (This potential 
coincides with the i/-component ($ = <& y ) of the complex self-dual electromagnetic vector 
potential = {0, 0, $ y , $ 2 }.) 



The symmetry reduced Einstein— Maxwell equations. It is well known that, for 
these space-times, the nontrivial part of the Einstein-Maxwell equations decouples into 
two parts. One of these contains the constraint equations which can be considered as 
determining the conformal factor f(u,v) (up to a multiplicative constant) in terms of the 



other functions: 



fu 

J 

fv 
f 



a. 



a v 



H u _a_ 
H 2a u 



H v a 



-ir^ |2 



\£u ~t~ 2 < D < D W 
|&, + 2$$, 



H 2 



(7) 



where all suffices mean the derivatives, and the right hand sides are expressed in terms of 
the complex Ernst potentials £ (u, v) and $(u, v) which are determined by the relations: 

Re£= ( Re$= A, 



d u (lm£) 
d v (lm£) 



a: 



l H 2 d u n + i($d u <f>-$d u $) 



<9 u (Im$) 
<%(Im$) 



-oC l H{d u A z - Qd u A y ) 
a~ 1 H{d v A z - ttd v A y ) 
(8) 

where A y , A z are the nonzero components of a real electromagnetic vector potential. 

The remaining part of the reduced electrovacuum Einstein-Maxwell equations are the 
dynamical equations for the metric functions a(u, v), H{u,v) and Q(u, v). It is very 
convenient to present these in the form of the hyperbolic Ernst equations for the Ernst 
potentials: 

(Re £ + $1) (2£ uv + ^£ v + ^£ u ^- (£ u + 2$$ u ) £ v - (£ v + 2$$„) £ u = 
(Re£ + $$) (2$ uv + ^ v + - (£ u + 2$$ u ) - (£ v + 2$$„) $ u = ( 9 ) 







These dynamical equations play the primary role in studies of colliding plane waves 
because they govern the nonlinear processes of the interaction between the waves. They 
also provide the integrability conditions for the constraint equations (|IJ) and the compat- 
ibility of the relations (JBJ) which relate all the metric functions and other field variables 
with the solutions of the Ernst equations ©. The constraint equations determine the 
remaining metric component - the conformal factor f(u,v). 



Matching conditions at the wavefronts. We say that a plane wave with a distinct 
wavefront, say u = or v = 0, propagates through a given background if any wave-like 
solution of © defined for u > or v > respectively matches appropriately on the 
corresponding wavefront to the given background solution which is defined in the region 
u < or v < respectively. 

For colliding plane waves, it is well known that when considering the junction con- 
ditions across the null hypersurfaces which correspond to the wavefronts, the familiar 
Lichnerowicz conditions, which require that there should exist a coordinate system in 
which the components of the metric and electromagnetic potential are at least of class C 1 
on the wavefront, should be relaxed in favour of the O'Brien-Synge conditions. These ad- 
mit also some other types of waves, such as impulsive gravitational waves or gravitational 
and electromagnetic shock waves. For the class of metrics defined above, these conditions 
imply that across the wavefront 

(i) the function a is of the class C 1 

(ii) the Ernst potentials £ and $ are continuous 

(iii) the conformal factor / is continuous 

It is important to note here (see below for more details) that the condition (iii), through 
the relations (|7j), may impose some further restrictions (additional to the conditions (i) 
and (ii)) on the behaviour of function a and the Ernst potentials near the wavefront. 
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3 Plane waves in a Minkowski background 



The Minkowski background. Within the class of metrics (JT|), the Minkowski space- 
time can be represented in the form 

ds 2 = dudv — dy 2 — dz 2 (10) 

or, in the notations introduced above, by the values 

a = l, £ = -1, $ = 0, f = l (11) 

(and therefore, H = 1 and Q = 0). Everywhere below, we shall consider this space-time 
as the background for the waves in the sense that was formulated in the previous section. 

Travelling waves. It is well known that the Ernst equations (jHJ) admit a wide class of 
plane wave solutions which can also be called "travelling waves" . For these solutions all 
potentials depend on one null coordinate u or v only and, therefore, correspond to plane 
waves travelling without any evolution in the positive or negative x-direction respectively. 
To distinguish these similar solutions, we shall call the waves depending on u, and trav- 
elling therefore in the positive x-direction, the "left waves" , and the waves depending on 
v and travelling in the negative x-direction the "right waves" . 

The left waves. For these waves all metric components and potentials depend only on 
u. It is very convenient to use the coordinate freedom u — > u(u) to specify the conformal 
factor f(u) = 1. Geometrically this means that the coordinate u is chosen to be an affine 
parameter on the null geodesies v = const, y = const, z = const which cross the wave. 
This class of waves can be described completely by the Ernst potentials 

{£(u), *(«)} (12) 

which are functions of the affine parameter u and which should satisfy only the signature 
condition 

H{u) = -Re£(u) - > 0. (13) 

The corresponding function a(u) is determined for u > as the solution of a linear 
differential equation which follows from the first of the constraint equations (J7|) for the 

case f{u) = 1: 



Hu , 1 
ti z 



|£ U + 2$$J 2 i 
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a = 0. (14) 



H 2 H 

with the initial data at u = which follows from the matching conditions. 

These waves for u > are combined with a Minkowski background for u < as 
illustrated in Fig.(l). These two regions are joined on the null hypersurface (wavefront) 
u = where all metric functions and the Ernst potentials (|12|) defined for u > (the 
wave interior region) should be matched with their Minkowski values 

The matching conditions (i)-(iii) of the previous section imply the following conditions 
on the wavefront u = 

£(0) = -l, $(0) = 0, a(0) = 1, a'(0) = (15) 



s 




Figure 1: The left wave propagating on the Minkowski background. The coordinate u is chosen 
to be an affine parameter on the geodesies v = const. In this case, / = 1 everywhere. 

(together with the conditions H(0) = 1 and f2(0) = 0). The last two conditions in (|15J) 
provide initial conditions for the equation (fb"f|) for a(u) which is then determined uniquely. 
It follows from the structure of equation (fT^j) for H(u) > that, for any nonconstant data 
£(u), the corresponding function a(u) decreases monotonically from its initial value 

a(0) = 1 with a'(0) = 0, when u runs from to some finite critical value u* where a(u*) 
vanishes. Since a is a measure of the area of two-dimensional sections of the null geodesic 
tubes, its decrease is a manifestation of the well known general focusing property of 
gravitational fields. The surface u = corresponds to the presence of a caustic, and 
represents a boundary of the space-time being considered. 

The right waves. Solutions in which the metric components and the Ernst potentials 
are functions of the null coordinate v only are exactly equivalent to those described above, 
except that they propagate in the opposite, negative x-direction. In an exactly equivalent 
way to that described above, we impose the coordinate condition f(v) — 1 so that v is an 
affine parameter and any solution is then described by the corresponding Ernst potentials 
{£(v), $(f)} and the function a{v) > satisfies an equation equivalent to (JHJ). Again, 
we introduce a wave of this type for v > by matching it with a Minkowski region for 
v < exactly as above. 

Other coordinates. Sometimes the wave solution simplifies considerably if we relax 
our choice of the null coordinate u (or v) as the affine parameter on the null geodesies 
crossing the travelling wave interior region. 

In practice, it is often convenient to avoid having to solve the equation (|14j) for a. 
Instead, the freedom u —>■ u{u) can be used to specify the form of a(u). However, this 
introduces a non-constant conformal function f(u) which, for any wave functions £(u), 
$(ii), can be calculated in quadratures from the constraint equations (J7|). Of course, the 
choice of the function a is restricted by the condition that it and its first derivative must 
be continuous on the wavefront and monotonically decreasing through the wave region. 
In addition, the choice of the Ernst potentials £{u) and $(w) should be also restricted by 
the condition for the regularity of the space-time geometry and the electromagnetic field 
components near the wavefront, and this imposes further restrictions on the character 



of a at the wavefront. This alternative approach will not generally be adopted in the 
remainder of this paper. 



4 Characteristic initial value problem for colliding 
plane waves 

Now we consider the situation in which two shock waves exist in the same space-time, 
are initially separated by some Minkowski region, and approach each other from oppo- 
site spatial directions. The left wave comes from the negative x-direction and has the 
wavefront u = 0, while the right wave comes from the positive x-direction and has the 
wavefront v — 0. 

Before their collision, their interior regions (u > 0, v < for the left wave and v > 0, 
u < for the right wave) are described completely by the Ernst potentials expressed as 
functions of the affine parameters on the null geodesies crossing these regions. These are 
represented by the functions £ + (u), for u > and S-(v), $_(t>) for v > for the 

left and right waves respectively. (For convenience, we denote functions associated with 
the left and right waves using the subscripts + and — respectively.) As discussed in the 
previous section, the function a which appears in the dynamical equations as an additional 
dynamical variable, is not an independent function. Rather, in the approaching waves, 
the functions a+(u) and a_(i>) are determined uniquely as solutions of the corresponding 
constraint equations with initial data determined by the matching conditions on their 
wavefront s. 

A collision of these waves at the point u = 0, v = gives rise to the existence in space- 
time of a wave interaction region (u > 0,v > 0). As in the approaching wave regions, 
this region also possesses a boundary which is determined by the condition a(u, v) = 0. 
However, unlike the boundaries u = and v = v* of the approaching wave regions, this 
boundary is generally a curvature singularity (see [3] for more details). The complete 
space-time is illustrated in Fig. (3). 

The problem of the construction of the solution in the wave interaction region for a 
given pair of approaching travelling plane waves is clearly a characteristic initial value 
problem. The matching conditions applied on the boundaries (u = 0,v > 0) and (v = 
0,u > 0) of the wave interaction region imply the continuity of the function a together 
with its first derivatives and continuity of the Ernst potentials on these characteristics. 
These conditions provide all the necessary characteristic initial values for these functions. 

First of all, it is easy to see that the characteristic initial value problem for the function 
a(u, v) separates from that for the Ernst potentials and its solution is trivial: 



because the matching conditions for the approaching waves with the Minkowski back- 
ground have been satisfied, and the functions a+(u) and a~(v) possess the properties 



In particular, at the point of collision we have the "normalization" a(Q, 0) = 1. 

We also have the characteristic initial data for the Ernst potentials in the interaction 




(16) 



a+(0) = a_(0) = 1 



a' + (0) = a'_(0) = 0. 



(17) 



region: 



£{u, 0) = e+(u) : 
$(w,0) = $+(«) 



8{0,v) = e-{v), 
$(0,«) = $_(«). 




in 




Figure 2: The collision of travelling plane waves propagating into a Minkowski background. The 
coordinates u and v in the wave regions before their collision are chosen to be affine parameters 
along the geodesies v = const and u = const respectively. The conformal factors in these regions 
are f+(u) = 1 and /-(f) = 1 respectively. However, in the wave interaction region, f(u,v) ^ 1 
and a(u, v) = a+(u) + a.-(v) — 1. 

These are identical to the functions which determine the structure of the approaching 
waves. The compatibility of these data at the point of collision follows from a correct 
matching of the original travelling waves with the Minkowski background. In this case, 
these data satisfy the conditions 

£+(0) = £_(0) = -1, $+(0) = $-(0) = 0. 

Thus, these also provide normalization conditions for the Ernst potentials at the point of 
collision: £(0,0) = -1 and $(0,0) = 0. 

Geometrically defined null coordinates (u,v) and (£,??). In the above formulation 
of the characteristic initial value problem, it seems appropriate to emphasize that two 
different pairs of geometrically defined null coordinates can be used. The coordinates 
(u, v) have been adapted to the orbit space of the isometry group where the dynamics of 
waves effectively takes place. These coordinates have been introduced in the approaching 
wave interior regions, prior to their collision, as the affine parameters on the plane wave 
rays and on the geodesies crossing the corresponding wavefronts. In the wave interaction 
region, these coordinates can then be defined uniquely using their continuity on the initial 
characteristics. However, in this region, these coordinates lose their property as affine 
parameters. 

It is also convenient to introduce another system of null coordinates (£,??). These 
are also defined geometrically but, unlike the (u,v) coordinates which are adapted to the 
orbit space of the isometry group, the definition of these coordinates is adapted to the 
internal geometry of the orbits themselves and, therefore, to the internal geometry of the 
wave surfaces. More importantly for us, the coordinates (£, rj) are better adapted to the 
analysis of the field equations. 

The relation between these two coordinate systems can be found using the definition 
(JH) and the expression (fTrjj) derived above. We can take the function j3, which is conjugate 
to a, as 

(3 = a + {u) — a~(v). 
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A constant of integration in this expression has been chosen to give /3 — at the point of 
collision. With this expression, £ and r\ have the forms: 

i(u) = 2a+(u) - 1, r)(v) = 1 -2a_(v). (19) 

Since «+(«) and a_(t> ) are monotonically decreasing functions, £ and 77 can be adopted as 
coordinates throughout the interaction region, with £ = 1 and 77 = — 1 on the wavefronts 
it = and t> = respectively. 

Unfortunately, when formulated in the coordinates (£,77), the characteristic initial 
value problem under consideration possesses some unpleasant properties. The problem 
arises from the fact that these coordinates, when extended to the initial regions, take 
the forms £ = 1 for u < and 77 = — 1 for v < 0, so that the expressions for the inverse 
coordinate transformations «(£) and v(rj) have nonanalytic behaviour near the wavefronts. 
This means that, even for physically acceptable solutions which possess no singularities on 
the wavefronts and which, therefore, are described by regular functions of the coordinates 
u and v, the metric components and potentials possess some singular behaviour near the 
wavefronts when expressed in terms of the coordinates £ and rj. In particular, although 
they are continuous on the wavefronts, they possess unbounded £- or ^-derivatives there. 

5 Monodromy transform and the integral evolution 
equations 

In our present construction, we use the general approach developed in [IS] . [T7] and [THj 
for the analysis of the structure and construction of solutions for integrable hyperbolic 
reductions of the vacuum Einstein and electrovacuum Einstein-Maxwell field equations. 
In this approach, which is called the monodromy transform (because of its analogy with 
the well known inverse scattering transform), every solution of the reduced field equations 
is characterized by two (for vacuum) or by four (for electrovacuum) functions which 
depend on the spectral parameter only. These are interpreted as the monodromy data of 
the corresponding fundamental solution of the associated linear system. The problem of 
the solution of the nonlinear field equations has been substituted in this approach by a 
solution of some linear singular integral equations whose scalar kernel is built from these 
monodromy data. 

A further development of this approach, which provides some new mathematical tools 
for the solution of characteristic initial value problems, has been presented in |2] - In 
that paper, the general solution of the associated spectral problem was represented by 
two scattering matrices. These dress the initial values of the fundamental solution of the 
associated linear system given on two characteristics which pass through some arbitrarily 
chosen point where the solutions are normalized. The monodromy structure of these 
scattering matrices is characterised by "dynamical" monodromy data whose evolution is 
determined completely by the initial data for the fields on these characteristics. In this 
context, the problem of the solution of the nonlinear field equations reduces to a solution 
of some linear integral "evolution" equations of quasi-Fredholm type. The scalar kernels 
of these equations are built from the dynamical monodromy data and the initial values of 
the fundamental solution of the associated linear system on the above-mentioned initial 
characteristics. In terms of solutions of these linear integral evolution equations, the 
corresponding solution of the nonlinear field equations can be calculated in quadratures. 
Thus, many properties which are revealed by these constructions are in principle well 



1? 



adapted for the solution of the characteristic initial value problem for various integrable 
reductions of Einstein's equations. 

As formulated in the cited papers, however, these constructions cannot be applied 
to a solution of the colliding plane wave problem. This is because the local analyticity 
of the solution is violated on the wavefronts. Moreover, as mentioned at the end of the 
previous section, although the components of physically acceptable solutions are regular 
functions of the (u, v) coordinates, their £- and //-derivatives have singular behaviour on 
the wavefronts. Fortunately, however, it has been found possible to derive a generalization 
of the previous construction which overcomes this difficulty This was reported briefly in 
PQ where the corresponding generalized linear integral evolution equations were presented. 
It was shown how the character of the singularities of the kernels of the quasi- Fredholm 
integral equations can be explicitly related to that of the characteristic initial data. This 
construction, whose detailed description and some applications we present below, remains 
valid and, eventually, opens the way for the consideration of a wide range of types of 
colliding waves with nonanalytical behaviour on their wavefronts. 



5.1 Associated linear system for the Ernst equations 

Consider a linear system for a complex 3x3 matrix function *&(£,r),w) which depends 
on two real null coordinates £, rj (as mentioned above) and a free complex (spectral) 
parameter w. To represent the electrovacuum Ernst equations, this system should be 
supplemented by additional conditions of the following two kinds. The first consists of 
algebraic restrictions on the matrix coefficients of the system. Together with the linear 
system itself, this can be presented in the form 



d £ V = U ^'^ N ^, rankU = l, trU 

9 r .* = V ^ ,T] \ y, rankV =1, tr V 
2i[w — rj) 



(20) 



The other conditions imply the existence of a Hermitian matrix integral of (|2(J|) which 
possesses the following structure 

Wj(iu) = WoM, 




4iO, n = -1 (21) 



where the Hermitian conjugation t is defined as ~W Q (w) = ~Wq(w), and W (u>) is an 
arbitrary (nondegenerate) Hermitian 3x3 matrix function of w. 



Equivalence of the matrix problem to the field equations. The conditions (|20p. 
(J21J) are equivalent to the symmetry reduced Einstein-Maxwell equations fT7\ I18j. For 
pure vacuum gravitational fields, the same set of conditions arise but for 2x2 matrices 
for which the third rows and columns of all matrices in (|2(J|) and (J21)) are omitted. In the 
fully general case, these conditions imply that the matrices U and V have the specific 
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forms 




-d v £, <%$), 



^4 



1 N 

P± i o 

v?± 1 



(22) 



in which p± = Cl ± — and g-j- 

Z2 



2<J> — 24> ( Q ± — ) . Also, the matrix W must he a. lineal- 



function of the spectral parameter w with the structure 

„2 



W = Ai(w - /3)fi + G, G 



^(fo 2 + 21; 

-4#fi- 
y -2$ 



) + 4$$ -4iffi - 4$$ 

4$$ 4# + 4$$ 

2$ 



-2$ 

2$ 
1 



(23) 



If in fl22j) - (|23|) . a, j3, H, Q, £, $ and $, are identified respectively with the previously 
defined functions a = |(£ — r]), (3 = |(£ + rj), the metric functions H and Q in ©, 
the Ernst potentials £ and $ = $ y , and another component $ = $ z of the complex 
electromagnetic potential, then (J2(J|) - (J21j) imply that all the field equations and relations 
given above are satisfied for these functions. The inverse statement is also true, so that all 
of the functions in ([22 )1 — (|23 |) can be calculated for any solution of the reduced Einstein- 
Maxwell equations. Thus, the set of conditions (J2*m) -(|21| ) for the four complex matrix 
functions 

¥(Z,7i, w), U(£,7?), V{Z, V ), W(£,7], W ) 

constitute a problem whose solution is equivalent to the solution of the symmetry reduced 
Einstein-Maxwell equations. However, there is a significant advantage in substituting 
an analysis of the matrix problem (J2(J)) - (J21|) in place of studying the nonlinear partial 
differential equations which govern the collision of gravitational or gravitational and elec- 
tromagnetic waves. It was just this step that was made in previous formulations of the 
monodromy transform approach, and we develop this further here. 



5.2 Normalization conditions 

Before we analyze the matrix problem (J2(J)) -(|21 )) . it is important to remove the gauge 
freedom which exists in these equations. Without loss of generality, we can impose nor- 
malization conditions for the values at some "initial" or "reference" space-time point of 
the matrix functions ¥ and W, as well as of the components of the metric and elec- 
tromagnetic potentials and the Ernst potentials. For the problem considered here, it is 
convenient to choose this point of normalization as the point at which the waves collide, 
i.e. with the point u = 0, v — 0, or £ = 1, 7] — —1. First, we recall that we have already 
normalized the metric functions and the Ernst potentials at u = 0, v = as follows: 

a(0, 0) = 1, H(0, 0) = 1, £{0, 0) = -1, 5(0, 0) = 
/5(0, 0) = 0, fi(0, 0) = 0, $(0, 0) = 0. 
For the matrix function St, we can adopt the normalization 

*(1,-1, W )=I. (24) 
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by using the transformation 77, w) —>■ 77, to)\f r_1 (l, —1, w). 

From (J2*T|) and (J2BJ), it is now easy to show that the above normalization conditions 
imply that the matrix integral ~Wq(w) in (J2*T|) satisfies 

W (ui) = 4iu;0 + diag (4, 4, 1) . 

5.3 The analytical structure of VE* on the spectral plane 

Everywhere below, \£ denotes the fundamental solution of the linear system (|2())l. normal- 
ized at the point of collision as in (124(1 . The general analytical structure of this solution was 
described in detail in JE], ^7j, ^H]- However, in those formulations, \I/ was normalized 
at some arbitrarily chosen space-time point at which all components of the solution of the 
field equations are analytic functions of the coordinates £ and 77. However, in the present 
case we specifically choose the point of normalization to be the point (£,77) = (1, —1) at 
which the waves collide. It is obvious that this point is not a point at which the field 
components are locally analytic. 



-1 v Z 1 

Figure 3: Cuts in the spectral plane w 

As in the regular case, with the normalization (|24|). the structure of the associated 
linear system implies that the fundamental solution 77, iu) (and also its inverse) pos- 
sesses four branch points on the spectral plane w. The relative positions of these points 
and the cuts L± which join them are indicated in Fig. 3. These singular points include the 
two "variable" points w = £ and w = t] which correspond to the poles of the coefficients of 
the linear system (|20|). As in the regular case, these singularities for ^ are branch points 
of order —1/2. In addition, there are two fixed branch points of & at w = 1 and w = —1 
which arise from the normalization (}2"4"|) imposed at the point of the collision. However, 
the order of the branching of ^ at these two points may differ from 1/2 as occurs in the 
regular case. Moreover, the character of these branch points is not universal in this case. 
In fact, it becomes dependent on the solutions and differs from one solution to another. 

Fortunately, we can proceed with our analysis without needing to specify the order of 
branching of ^ at the points w — 1 and w = — 1 on the spectral plane. Following mainly 
the same simple arguments given in [THl HZ] and the ideas suggested for the regular case in 
j2] , we present below a generalization of this construction that is available for the analysis 
of fields which have nonanalytic behaviour on the characteristics. 

Because the rank of the matrix coefficients of the linear equations (|2~Uj) is everywhere 
not more than 1, near each of the singular points of only one of three linearly inde- 
pendent solutions of the associated linear system can be branching, while the two others 
should be holomorphic. This immediately implies that the normalized fundamental solu- 
tion ^ in general has the local structure on the cuts L± which can be expressed in the 
forms: _ 

L + : 77, w) = o|> +(£, 77, w) ® k+H + M+(f , 77, w) 

(25) 

L_ : T),w) = * _(£, v, w) ® k_(tu) + M_(£, 77, w) 

where, as in the regular case, the coordinate independent components of the row vectors 
k + (w) and k_(w) are regular near the cuts L + and L_ respectively or, at least, are 
not branching. However, they can have isolated singularities at the endpoints w — 1 

IF) 



and w — — 1 respectively. The same is true for the matrices M + and M_ on L + and 
L_ respectively. In the regular case, the vectors ip ± (£, 77, w) always have the structures 

$ + (£>V,w) = ^J^^+(C,V, W ) and ip ^,rj,w) = y^^tp _(£, 77, w) where the vectors 
ip + and ip _ are holomorphic on the cuts L + and L_ respectively. In the general case, the 
components of the column vectors \p ± (£, rj, w) also have branch points at the endpoints 
of the corresponding cuts L±. However, the characters of their singularities at the points 
w = 1 and w = —1 respectively can be different for different solutions. Their specific 
characters are determined by the physical properties of the wavefronts as incorporated in 
the characteristic initial data for the approaching waves. 

5.4 Characteristic initial value problem for VE* and analytical 
structure of its initial data 

Following a similar construction as for the regular case given in 0, we now reformulate 
the characteristic initial value problem for the hyperbolic Ernst equations described in 
section 0] as an equivalent characteristic initial value problem for the matrix function St. 
Later it will be shown that the solution of this problem can be reduced to a solution of 
some simple linear integral equations. 

First, it is easy to observe that the given initial data ([18)1 on the characteristics en- 
able us to calculate the matrix coefficients in (j20j) on these characteristics, i.e. we can 
determine the values of U(£, 77 = —1) and V(£ = 1,T]). This enables us to introduce two 
matrix functions w) and \&-(j],w) which are the normalized fundamental solutions 

of the associated linear system on the characteristics £ = 1 and r] — — 1: 

I *+(i, W )=i i *_(-i,«,)=i 

(26) 

The coefficients of each of these systems are completely determined by the initial data for 
the fields on the corresponding characteristic. Thus these matrices form the characteristic 
initial data for the required solution for \I>(£, 77, w): 

w) = *(f , -1, w), *_(r/, w) = r/, w) 

Applying now the same arguments that were used in El > we conclude that the 
boundary values w) and ^-(77, w) should possess analytical structures on the spec- 

tral plane which are very similar to that of T],w). Namely, *$?±(w = 00) = I and 
ty + and its inverse are holomorphic outside L + , while ^_ and its inverse are holomorphic 
outside Moreover, their local structures on the cuts are given by the expressions 

L + : * + (^w)=$ 0+ (^w)(g)k + (w) + M 0+ (^w), 

(27) 

L_ : (77, w) = ip _ (77, w) <g> k_ + M _ (7/, to) 

where h±(w) are the same as for ^ in (|2~5|l . and Mo+(^, w) and Mo_(?7, w) are non- 
branching on L + and L_ respectively. However, ty± may have isolated singularities at 
the endpoints w — 1 and w = —1 of these cuts respectively. As for the general case 
described above, the vectors tp 0+ and \p _ have branch points at the ends of the cuts L + 
or L_ respectively. Also, the character of their branching at these points is completely 
determined by the initial data. 



Ifi 



5.5 The "scattering" matrices and their structures 

We now introduce, in analogy with the regular case, the "evolution" or "scattering" 
matrices x±(£,iVi w )i which represent T],w) in two alternative forms 

V(€,7i,w)=x+(€,V,w)--&+(£,w), ^(^rj,w)= X -(^V,w)-^-(rj,w). (28) 

To understand the analytical structures of these matrices on the spectral plane, we express 
them using the above relations in the form 

x +(e,^,w) = *(e,r?,w)-*; 1 (e,w), x-fc^w)^*^™)-*: 1 ^,™)- (29) 

It may be noted that these matrices are solutions of linear equations with well defined 
initial conditions. We obtain these equations for x ± after multiplication of the second of 
two linear systems in ([20|) by \I>7_ (£, to) and the first one by VfC (77, to) from the right. 
We then find that each of the matrices x + and x - are determined by the linear systems 
with obvious initial conditions 



2t(w-0 A ~ ^ ^ + 2i(w- V ) " + (30) 
X-(l,r],w) = l ( x+(C,-l,w) =1 

For any given local solution of the Ernst equations each of these two systems admits 
a unique solution. In both of these cases, the structures of the linear systems and the 
initial conditions show immediately that the corresponding fundamental solutions x+ or 
X- should be holomorphic on the spectral plane w everywhere outside the cut L_ or L + 
respectively. Each of these fundamental solutions possesses only two branch points, which 
are located at the endpoints of the corresponding cut, and a finite jump on this cut. In 
addition, it is easy to show that these matrix functions should also possess the property 
X±(£,t?,w = 00) = I 

All of the above features enable us to represent each of these matrix functions as the 
Cauchy integrals of their jumps over the corresponding cut: 

x+ ^ v ,w) = I + - [ -M^-dC- x-^V,w) = I+- [ -M^dC + (31) 
m J — u> in J ( + — w 

L_ L + 

where [. . .]^ is the jump (a half of the difference between the left and right limits) of a 
function at the point w = ( on a cut. It may be noted that we usually choose directions 
on the paths L + and L_ shown in Fig. 3 from w = 1 to w = £ and from w = — 1 
to w = 77 respectively. However, the above definition of the jumps makes the integrals 
()31j) insensitive to the choice of these directions. It is also necessary to note here that 
the convergence of the Cauchy integrals in (|3T|) at w = —1 and w = 1 respectively has 
been assumed. This conjecture has been proved for the regular case, and is confirmed 
for all known examples which possess different degrees of non-smooth behaviour on the 
characteristics. 

Let us now explicitly calculate the structure of the jumps of x± on the cuts using the 
definitions the local structures (|23j) as well as the analyticity of and its inverse 
on L_ and and its inverse on L + . Simple calculations lead to the expressions 

[x+] r _ = [tj>_] T - ®m_(£,r_), [x-] r+ = [tp+]r+ ®m + (r?,r_), (32) 
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where r + G L + and t_ G L_. Substituting these into the integrands of ([31)1. we obtain 
x+U,^,w) = I+— / d(- 



ni J — w 

,t , Tj _ 1 /- [tP + ] C+ ^m + faC + ) (33) 

These expressions involve a new, evolving kind of vector-functions 

m + (?7, w) = k + (u>) • ^I 1 ^, w), m_(£,w) = k_(w) • (34) 

These "dynamical" monodromy data replace the coordinate independent vectors \s.±(w) 
which arose earlier in similar integral representations for St normalized at a regular point. 
(The interpretation of these vector functions was explained in [2] for the regular case.) 

5.6 Integral "evolution" equations 

The alternative representations of \I/ that are given in (|28J) implies an obvious consistency 
condition 

X+(f,»7, «0 = X-fov,w)*4v,w)- (35) 

It is useful to recall that the matrix functions and x - are holomorphic on the spectral 
plane everywhere outside L + , while ^_ and x+ are holomorphic outside L_. In addition, 
the values of these four matrices at w = oo are all equal to the unit matrix I. Therefore, 
the right and left hand sides of ()35j) are analytical functions outside the cuts L + and L_. 
Hence, for the condition (|35|) to be satisfied, it is necessary that (here and everywhere 
below, r + G L + and r_ G L_) 

X+(t + )[9 + ] t+ = [ X -] r+ *-(r + ), |x + ] T _* + (r_) =x-(r_)[*_] T _. (36) 

Here (and in some expressions below) we omit for simplicity the parametric dependence 
of the functions on the coordinates £ , r], or both £ and r\. 

We now substitute into (|3Tjl the integral representations (JSHj) . and the expressions 
The expressions for the jumps of on the cuts L± can then be derived easily from the 
local representations (|2*7|) : 

[*+]r + = [^o+]r + ®k+(r+), r+GL+ [*_] T _ = [ip _] T _ ®k_(r_). (37) 

It is convenient to denote the jumps of the column-vector functions \\> ± on the cuts L± 
as 

<t>+(€,V, T +) = [*+]t+, c|>_(£,J7,t_) = [xF_] T _, 
and the initial values of these jumps as respectively 

<t>0+(£,^+) = [^0+]r+, <|>o_(77,r_) = ft> ~]r_- 

Substituting the expressions (J3~3*)) . (J32*|) and (}3T|) into (JSTIj) . we arrive at the following 
coupled pair of linear integral equations for the functions rj, r+) and t|> _(£,??, r_), 

which are defined on the cuts L + and L_ respectively: 

* + (r+) - J S + (r + , C_)d> _(C-) dC- = * o+(r+) 
|5_(r_,C + )c t)+ (C + )rfC + = *o-(r-) 



(38) 



*_(T_)- 
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where t + ,(+ £ L + and r_,C_ G L_. Also the scalar kernels of these equations take the 
forms 

fm_(£,C-) • *o+(£,T + )J (m+(??,C+) • *o-(r/,r_; 

T+ ' C - } = i7 r(C_ - r + ) ' ^ T -' C+) = ,vr(C + - r_) " 

(39) 

The equations ffTTH^) can easily be decoupled into two separate equations for the vec- 
tor functions cf> 77, r + ) and <t> _(£, 7/, r_). Substituting c|) _(£, 77, r_) from the second 
equation in (J3~%|) into the first, and inversely, we obtain the decoupled integral equations 



- J F + (T + ,( + )4> + (C+)d( + = f+(r 4 



(40) 



in which the kernels and right hand sides are given by 



(41) 



L- 

F-( r -,C-) = J 5-(r_,x+)5+(x+,C-)^X+, 

f+(r + ) = * 0+ (r + ) + y 5+(r + ,X-)*o-(X-)dx-, 

f_(r_) = *o-(r_) + y 5_(r_,x+)c()o+(x+)rfx+- 

However, although they are decoupled, these equations possess more complicated struc- 
tures than those in (|3~H|). 

All the coefficients of the integral equations (f3~H|). and the equivalent decoupled pair 
are determined by the initial data functions w) and ^f^(r],w). These are 

determined, in turn, by the characteristic initial data for the field variables (e.g., the 
Ernst potentials). The relations of the coefficients with the characteristic initial data will 
be clarified in more detail in subsequent sections where we consider particular examples 
of the calculation procedure. 



6 Principal algorithm for solution of the problem 

6.1 Specification of the initial data 

As explained in sections El and EJ we may consider the initial data to consist of the 
Ernst potentials, given as functions of the affine parameters u and v along the initial null 
characteristics 

{£+(«),*+(«)}, {£-(«),*-(«)}. ( 42 ) 
Alternatively, the initial data may consist of the metric functions and the complex elec- 
tromagnetic potential (again as functions of the affine parameters u and v) 

{H + (u), n+(u), $+(«)}, {#_(«), $-(v)} (43) 



1Q 



In either case, we require the continuity conditions 



£t(0) = -l, $±(0) = 0, H ± (0) = 1, n ± (o) = o. 



(44) 



Because of the choice of u and v as affine parameters on the initial characteristics, for the 
conformal factor /, we have f+(u) = 1 and f-(v) = 1. Also, the function a(u, v) is given 
explicitly by pfijl . where the functions a + (u) and a_(t>) are determined from the linear 
equations with the well defined initial conditions: 



Q_ 



H' a 



a+(0) 



1. 



<(0) 



#4 



+ 41$ 



/ |2 



(45) 



0. 



where H + {u) = — Re£ + (w) — $ + (m)$ + (m) and 



a" - «' 



+ 



a_(0) = l, a'_(0) = 



£> +2$_$' I 2 , ,„ 

^-+4$'_ 2 



(46) 



where H_(v) 



— Re£_(v) — $_(u)$_(u). Having obtained the functions a + (w) and 
a_(u), the functions and 77(f) must now be determined using the definitions (|4"|). 

A nice feature of the above specification is that the functions f!42j) or (jjHj) can be taken 
as initial data functions and chosen arbitrarily, provided only that the conditions f)44j) 
are satisfied. In this case, we have to solve the linear equations (}43j) and (j4T?j) to find the 
corresponding functions a + {u) and a_(w). 



6.2 Characteristic initial values for ^ 

Everything is now ready for the calculation of the characteristic initial values for *&(u, v, w). 
From the initial data, we must calculate the matrix coefficients U(£(tt), — 1) = U(u, v = 0) 
and V(l, 77(f)) = V(ii = 0,v) of the linear systems (j^j) . using the definitions (j^jl and 
expressing them in terms of the u, v coordinates. We should then solve the equations ()26)1 
and determine the normalized fundamental solutions & + (u,w) and *&-(v, w). 

It is clear, of course, that this step is crucial for a practical realization of the method. 
In practice, however, for a great majority of the initial data, it is found to be impossible 
to solve these equations explicitly in a closed form. Nonetheless, it may be argued that 
there should exist some infinite hierarchies of the initial data for which the solutions of 
these linear systems can be found in a desired form. In this case, it may be hoped that 
such cases correspond to sets of initial data which are reasonably dense in the whole space 
of the initial data and at least contain cases that are of particular physical significance. 



6.3 Construction of the integral evolution equations and their 
solution 

Once the matrix functions ^f + (u,w) and ^f-(v,w) have been found explicitly, we are 
able to calculate the fragments of their local structures (|27jl . i.e. the monodromy data 
vectors \c±(w), the dynamical monodromy data vectors m + (v, w) and m_(u, w) as well as 
the characteristic initial values ip 0+ (u, w), ip q-(v, w) of the functions tjj + (u, w), ip _(v, w) 



and the initial values of their jumps <p o+(u, w = r+), o-(v, w = r__) on the cuts L + and 
L_ respectively. This is all we need to construct the kernels and the right hand sides 
of the main integral evolution equations (J3%|) or the kernels and right hand sides of the 
decoupled form (jjUj) of these equations. 

It is then necessary to solve the linear integral equations (J35j) or (J4*U|) with the corre- 
sponding kernels and right hand sides. This is the second crucial step in our algorithm 
because, even in the case in which we are able to find *& + (u, w) and *S? + (u, w) explicitly, 
there are no guarantees that the solution of the linear integral equations (fHH|) or (|4()jl can 
be found in an explicit form. However, if the particular forms of the coefficients of these 
integral equations (derived for the specified initial data) enables us to find their solution 
explicitly, then the Ernst potentials and all the components of the metric and the electro- 
magnetic vector potential in the wave interaction region can be expressed in quadratures 
(or even explicitly) in terms of these solutions. 



6.4 Calculation of the field components 

As in the regular case j2j, all components of the solution of the characteristic initial value 
problem can be expressed in terms of the matrices R±(u, v) defined by the asymptotic 
expansions of the scattering matrices 

X±(u,v,w) = I + -R ± (u,v) + 0(—). (47) 
w w z 

The integral representations suggest the following expressions for the matrices R±(w, v) 
in terms of solutions of the integral evolution equations ()38|) or (|40p and the dynamical 
monodromy data vectors m + (t>, r+) and m_(w, r_) 

R + {u,v) = — - / $_(u,v,C-) ® m_(w,C-)dC- 
iri J 

1 L ~r (48) 

R-{u,v) = ; / 4> + (u,v,(+) 

m J 

Using these expansions in (j20j) and (|2T|) leads to the expressions for the matrices U and 
V: 

V(u,v) = U(u,0) + 22<%R+, V(u,v) = V(0,u) + 2za„R_, (49) 

where we have to transform the £- and ^-derivatives into u- and f-derivatives respectively. 
In this way we also obtain the alternative forms of the matrix W 

W{u, v, w) = W{u, 0, w) - 4z(fJR+ + R^n) = W(0, v, w) - 4i(OR_ + Rlfi) (50) 

where the constant matrix O was defined in (j21jl . In accordance with (|23jl . the values 
of the matrices W(w, 0, w) and W(0, v, w) in (jBTIj) are determined completely by the 
characteristic initial data for the fields. 

The expressions ()49)1 provide us with the following alternative forms of the expressions 
for the Ernst potentials: 

£ (u, v) = £(u, 0) — 2i(e 1 ■ R + (u, v) ■ e 2 ) = ^(0, v) — 2i(e x ■ R_(n, v) ■ e 2 ) 

(51) 

v ) = $(«, 0) + 2i(ei • R+(u, u) ■ e 3 ) = $(0, u) + 2i(ei ■ R_(u, u) • e 3 ) 
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where ei = {1, 0, 0}, e 2 = {0, 1, 0} and e 3 = {0, 0, 1}. It may be helpful here to mention 
also some additional interesting relations for the matrices TL±(u,v), such as 

£(u, 0) = -1 - 2i(ei • R_(m, 0) • e 2 ), 0) = 2i(ei ■ R_(it, 0) • e 3 ), 

(52) 

£(0, u) = -1 - 2i(ei ■ R+(0, v) ■ e 2 ), $(0, u) = 2i(ei • R+(0, u) ■ e 3 ). 

For later reference, we give here also the expressions (j51j) in a more explicit form 

£(u,v) = £+(u) + - [ (e x • (it, t),C-))(m_ (it, C-) ■ e 2 )dC- 



7T_ 



(53) 



^ (ei • c)j + (M,f,C+))(m + (?;,C+) • e 2 )dC+ 
<$>(u,v) = $+(«) - ^ / (ei • *_(«,v,C-))(m_(«,C-) ■ e 3 R- 

L- 

where we have returned to our basic notations, with £+(u), $+(it) and 8-{v), $_(f) 
denoting the characteristic initial data for the Ernst potentials on the characteristics. A 
more detailed description of the various steps of this algorithm is given in the following 
section, where we present various particular examples. 



7 Examples of solutions 

7.1 Collision of two gravitational impulses 

As a first example, let us consider the collision of two impulsive gravitational waves. 
This case is well known. It is described by the solution of Khan and Penrose ^H] if the 
polarizations of the approaching waves are aligned. Its generalization to the case when 
the waves possess nonaligned polarizations was found later by Nutku and Halil 20J . We 
consider these cases here to illustrate that our algorithm works in this situation without 
introducing any essential difference between the linear (aligned) or nonlinear (nonaligned) 
cases. 



Gravitational impulses before their collision. The gravitational waves under con- 
sideration have impulsive components located prior to their collision on two characteris- 
tics, and the space-time behind these is again part of a Minkowski space-time. Although 
the metric functions are continuous on the wavefronts, some of their derivatives in direc- 
tions which cross the wavefronts are discontinuous and the curvature possesses a distribu- 
tional character. In this case, the equation PB^l for the left wave as well as the equation 
()46|) for the right wave should be solved together with the conditions that all projections 
(JOJ) of the Weyl tensor should vanish behind the wavefront for each of the approaching 
waves. The solutions subject to the conditions f)44j) take the forms 

a + {u) = 1 — k 2 + u 2 a_{v) = 1 — k 2 _v 2 , 

£ + {u) = -l-2k + e~ i5 +u- k 2 + u 2 , £_{v) = -1 - 2k_e~ i5 -v - k 2 v 2 , (54) 

$»=0, /+(«) = !, $>)=0, /» = 1, 
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where k± and 5± are arbitrary real constants. The corresponding expressions for H±(u) 
and Q±(u) are 



HAu) = 1 + 2uk + cos 5, + k 2 u 2 



HAv) = 1 + 2wA;_ cos <5_ + ArV 



n + (u) 
njv) 



2uk + sin <5 + 
2vk_ sin 5_ 



(55) 



It may be noted that the constants k ± play the role of amplitudes for these waves. They 
determine the scales (distances) in the affine parameters u and v at which the waves can 
focus initially parallel null geodesic rays which cross them. The parameters 5± determine 
the polarization of the waves. It is clear that either 5 + or S_ can be transformed to 
zero by a rotation of the ignorable coordinates y and z. However, it is the difference 
5 + — <5_(mod27r) which has physical significance, determining the relative polarizations of 
the approaching waves. However, we keep below both the parameters 8± to maintain a 
symmetry of the subsequent expressions. 

Using (JT5J1 and ([54)1 . we obtain the functions a and /3 in the interaction region as 



a(u, v) — 1 — k 2 u 2 — kA) 2 



(3(u,v) 



2„.2 



k_v - k z + u 



From this we obtain the relation between the null coordinates £ and rj and u and v 
respectively: 

£ = l-2k 2 + u 2 , rj = -1 + 2k 2 v 2 . (56) 



Calculation of the "in-states" *Sf + (u, w) and &Av,w). With the characteristic ini- 
tial data (J54")) together with we can calculate the matrix coefficients U(£, r) = —1) 
and V(£ = 1, rj) of the linear systems (J26|) on the characteristics rj = —1 (i.e. v — 0) and 
£ = 1 (i.e. u = 0) respectively: 



U(u, 0) 
V(O^) 



•^K^-M)]®(i,<,o) 



2ik_v 



—i(e 



( "'- + k_v 
5 - - k_v) 




(57) 



-z\ 



Substituting differentiation with respect to u and v rather than with respect to £ and 
rj in (|2l)j) and using the coefficients (J57)l . these equations can be solved explicitly. The 
normalized fundamental solutions take the forms: 



A 



-i 



2(w- 1 



a: 1 



w - 1 - 2k + e~ i5 +u 
-i(w - 1 + 2k+e- i5 +u) 




ty + 1 + 2k_e 



. i (to + 1 - 2k_e~ i5 -v) 



1, z, 



1, -z, + 



/I 


z 


o\ 


2 


~2 


z 


1 











2 


2 


\o 





1/ 


/ I 


/ 


o\ 


2 


2 


i 


1 




~2 


2 





V o 





1/ 



(5? 



while the inverse matrices are 



—i 



(1 i 0) + 



i(w-l- 2k + e~ iS +u) \ 

2{w - 1) 
w-1- 2k 4 .e~ iS +u 



* 1 = -X. 
2 



1 

.0, 



(1 -i 0) 



/ w - 1 + 2k + e~ ld +u 

2(w - 1) 

i(w - 1 + 2k + e~ iS +u) 

2(w - 1) 2(w - 1 

V 1/ 

/ iu + 1 - 2k_e~ l5 -v i(w + 1 + 2k_e~ l5 -v) \ 

2(w + 1) 2(w + 1) ° 

i(w + 1 - 2k^e~ i5 -v) (w + 1 + 2k_e~ i5 ~v 



2(w + l) 




2(iu + l) 






V 1/ 

(59) 

where the functions \ + (u,w) and \_(v,w) are determined by the expressions (see also 
(33)) 



A + (w,w) 



w — 1 V w + 1 

which satisfy the additional conditions A + (u, w = oo) = 1 and \_(v,w = oo) = 1. 



(60) 



Construction of the integral evolution equations. For our construction, we now 
need to determine certain fragments of the structures of these matrix functions on the 
cuts L± defined in (127)1 . These are the conserved monodromy data vectors 

k + H = (i, 2,0), k_H = (i, -i,o) 

and the initial values for the vector functions whose jumps on the cuts L + and L_ respec- 
tively are the unknowns in the integral evolution equations: 



ip 0+ (w,w) 



A: 



2 (w-i: 



a: 1 



2(w+ 1 



w-1- 2k + e~ iS +u 
-i(w - 1 + 2k + e~ i5 +u) 


w + 1 + 2k_e~ i5 -v 
i(w + l- 2k_e~ i5 -v) 



We now calculate the dynamical monodromy data vectors using the definitions 
giving 

, , s A 2k_e- i5 -v 2k_e- i5 ~v, n ^ 

m + (v,{ + ) = ( 1- - — , ), 



m_(u, C_ 



1 + 



C+ + 

2fc + e" i,5 +M 



C+ + 

i(i- 2fc ; e :^), o 



C--1 ' v C- 

and the kernels S± defined in f!39j) . for which we obtain the expressions 



S+K V,T + ,C-) 



1 2fc + e- i<5 +u[A; 1 ] T , 



S , -(i*,v,t_,C+) 



1 2k_e~ i5 -v[\Z 1 ] T _ 



m ( C _ _ 1)(r+ _ i) > -v- - ^ OT ( C+ + 1)(r _ + i) 
For this case, the kernels JF ± and the right hand sides f± of the decoupled integral evolution 



equations (pfUj) denned in (JlTj) take the forms 

1 2k+k_e i( -- s +- s -^ 



f + (r + ) = ^IK'U 
where u = a/1 — k\u 2 and v = v/l — Pt> 2 . 



27T u( r+ -l)(C + +l) 



[ai 1 : 



1 2k+k_e i (- 5 +- 5 -1 



UV 



m w(r_ + 1)(C- - 1) 



1 

— z 

1 

/ 



2fc.e" 



- u 



[A: 1 ],. 



1 + k_e 



v(t+-1) \i{l-k_e~ ld -v) 



2k_e 



-i<5_ 



u(t_ + 1) - k + e 



-i5. 



u 
u) 



Solution of the integral evolution equations. The equations (jHHjl and (jjDj) show 
unambiguously that the unknown functions ± should possess the following structure: 



To. 



i 



i 



^4+ 

i \b + 



[A 



-ii 



if 1 

2 W 



1 



i V#- 



where A±. and -B-t are unknown functions of u and t>. Substituting these expressions 
into (|38|) and explicitly calculating the integrals using analytical continuations and the 
elementary theory of residues leads to algebraic equations whose solution is 



A 4 



fc + e~ +«(« + k_e v) 
uv — k + k_e l (~ s +~ s -^uv 
ik,e~ l5+ u(u — k_e~ l5 ~v 



uv 



- k + k_e i (- 5 +- 5 -1 



A. 



k_e v(v + k + e +u) 
uv — k + k_e { (~ s +~ s -)uv 
ik_e~ lS ~v(v — k + e~ %s+ u) 



uv 



uv 



- k+k_e i (- 5 +- 5 - s > 



uv 



Calculation of the Ernst potential. It remains now to substitute this solution into 
the general expression (J53j) for the Ernst potential S. Calculating the corresponding 
integral leads to the expression 



£(u, v) 



kW - k' z v 



2„.2 



2(u + k_ve~ iS -)(v + k + ue~ i5 +) 



uv 



- k,k_e^- 5 +- 5 -^ 



uv 



It is not difficult to check directly that this solution satisfies all of the characteristic initial 
conditions of the problem under consideration and that it coincides with the corresponding 
expression for the Nutku-Halil solution. The case 5 + = <5_ = can be seen to correspond 
to the Ernst potential of the Khan-Penrose solution. 

7.2 Collision of gravitational impulse with electromagnetic step- 
like wave of constant amplitude and polarization 

As a second example, let us now include a step electromagnetic wave with one of the 
impulsive gravitational waves as described in the previous example. Specifically, let us 
consider a step electromagnetic wave to come from the right such that, behind the wave- 
front, it has constant amplitude and polarization. This is expressed by the condition 



const, 



(61) 



where the nonzero value of 0o determines the constant amplitude and polarization of 
the electromagnetic wave, and the vanishing of \l/o indicates that there is no associated 
gravitational wave component behind the wavefront. 

It is first necessary to determine the initial metric functions which give rise to the 
specific components (ffTT]). These can be determined using the expressions (JHJ), together 
with the field equation (pftjj) and the definition of the Ernst potential ©• It is found that 
a general solution of these equations, which satisfies the junction conditions (J4lj) is given 
by 

a _(v) = l-(l + k 2 _) sin 2 (£_v), f_(v) = 1, 

g_(v) = -1 - 2k_ sin(£_w) - 2k 2 (l - cos(£_v)), 

$_(t>) = e^- [sin(£_v) + k_(l - cos(£_v))], (62) 
H_ (v) = [ cos(£_v) + k_ sin(£_v)] \ = 0, 

= -i e *7- [sm{£_v) - k_(l - cos{£_v))\ . 

In the above expressions, £_ and 7_ are real constants which determine the values of the 
constant amplitude and phase of the electromagnetic field such that the above scalar is 
0o = — 2£_e t ' y ~ . In addition to these parameters, a nonzero value of the real constant 
parameter k_ represents an impulsive gravitational wave localized on the wavefront. If we 
want to consider a purely electromagnetic wave, we simply put k_ = 0. However, in this 
subsection we do not use this simplification, and permit the parameter k_ to be nonzero. 

It is now necessary to obtain the geometrically defined coordinates in the interaction 
region. Using (|16|) with (|54|) and (}62|) as the characteristic initial data, we obtain 

a{u, v) = 1 - (1 + k 2 ) sin 2 (£_v) - k 2 + u 2 , /3{u, v) = (1 + k 2 ) sin 2 (£_w) - k 2 + u 2 . 

This leads to the following expressions for £ and r\: 

i = \-2k\u 2 , rj = -I + 2(1 + k 2 _) sin 2 (£_v). (63) 



Calculation of the "in-states" *& + (u,w) and ^^.(v,w). Using the characteristic 
initial data (|53j), (jB3|) and (jHH), we can again calculate directly, first the matrix coeffi- 
cients U(£,?7 = —1) and V(£ = 1,7/), and then the corresponding fundamental solutions 
*& + (u, w) and ^^(v,w) of the linear systems on the characteristics r\ = — 1 (i.e. 
v = 0) and £ = 1 (i.e. u = 0). However, since the left wave in this case is the same as 
that of the first example, the expressions for *ff + (u, w) and its inverse are given by ()58|) 
and (}59|) . The result of similar calculations for *&_(v , w) and its inverse for the right wave 
is 

ty_(v,w) = lp _(f,W 

1 I k_ sin(£_v 




M -(v,w) 



2ik_ -2k_ 0, 



and for the inverse matrix we obtain 

/ 2i 

= L(w) ® (p -(v,w) + N -(v,w), L(w) = (w + 1) ( 1, z, - 1) 



«P0-(v,w) = --, , „ w r ( 1, -i, 



(w + l)(w-wi-) V ' 2A; 
sin(£„t>) / ^" . 



k 2 cosUv) I 1 \ / z \ / 1 00 

\-4i*_y ^ 2A; -^ -2*- 0, 

where u>i_ = 1 + 2k 2 , and the function A_(w, to) was defined in (fHUJ) . 



Construction of the integral evolution equations. We now need the conserved 
monodromy data vectors 

k + H = (i, i, o) , k_H = (1, -i, -^j 

and the jumps <p 0+ and <$> _ of the column-vector functions ^ 0+ at r + G L + and \p _ at 
r_ G L_ respectively: 



*o+0,T + ) = [A + X ] r+ 



*o-(u,r_) = [A: 1 ]^ 



















| k + e-' l5 +u | 













sin(£_t>) I \ k 2 cos(£_v) 




t_ + 1 \ n / r - — w l- 









, : 




\2z'A; 



which are the characteristic initial values for the unknowns in our integral evolution 
equations. Then, we obtain the dynamical monodromy data vectors using (|H4)I and the 
above expressions for the matrix functions *& + (u,w) and &_(v,w) and their inverse: 

/ >- \ ( ^ ■ r\ 2k_ sin(£_f . . 
m + (v,( + ) = 1, i, . \ M l, 



C+ + 1 V"' ' 2k 
( i \ 2k,e~ lS+ u ( \ 

m.(«,c-) = (i,-i,- 5r )+-^ri-(i,*,o) 

Then, for the kernels S± defined in (|39|) . we obtain the expressions 

q / t , 2fc + e- tf + M [A; 1 ] T+ 2£;_sin(^)[A: 1 ] T _ 

S + (W,T + ,C-) = --— 777 7T, 5_(U,T_,C+)- 



<7r(C--l)(r + -l)' v ' zvr(C + + l)(r_ + l)- 

We now have to solve the integral equations (jHEJ) for the functions 4> + and (\> _ using the 
above expressions for the kernels and the characteristic initial values c)j 0+ and c£> _ . The 
procedure involved can be reduced to the solution of an algebraic system if we observe 
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that the structure of the kernels calculated above and the right hand sides of these inte- 
gral equations imply the following dependence of the unknown functions on the spectral 
parameter: 



* + = IK'U { a + + -5^-, } , * _ = [AlV | A. 



B 



+ 



c 



w. 



where the vector coefficients A±, B± and C_ are independent of r ± , although they can be 
functions of the coordinates u and v. Substituting these expressions back into the linear 
integral equations, we obtain that the coefficients A± and C_ should be expressed as 



A, 




A 




C = -k 2 cos(£_v) 




whereas the other coefficients should satisfy algebraic systems whose coefficients are ex- 
plicitly calculated integrals and whose solutions are 



B 



B 



u 



Z(u, v) 



k_ sin(£_v) 
Z(u, v) 



/ ucos(£_v) + k_ sin(£_v) \ 

iu cos(£_v) — ik_ sin(£_v) 

\ 4ik_u(l - cos(£_v)) ) 
/ k + e~ l5+ u cos(£_v ) + v \ 

ik + e~ l5+ u cos(£_f ) — iv 

\4ik + k_e- i5 +u(l - cos(£_v))/ 



where we have used the notations u = \/l — k\u 2 , v = y 1 — (1 + k 2 ) sin 2 (£_t>), and 

Z(u, v) = uv — k + k_e~ lS+ usin(£_v) . This solution of the integral evolution equations, in 
accordance with (jHBJ), leads to the following expressions for the Ernst potentials: 



£ (u, v ) = — 1 — k 2 + u 2 



2e- i5 +k + ucos(£_v) 



, n 2k_ sm(£_v) 
cos(£_v)) —, — s - (v + e 



Z(u, v) 



'k+u) 



Z(u, v) 



<£>(«, t>) = k_(l - cos(£_v)) + 



(u + k_ sin(£_v)) 



(v + e %5+ k + u cos(£_f )) sin(£_f ) 
Z(u, v) 



This is a new family of solutions which, for particular values of its parameters, reduces to 
various known solutions. Specifically, for pure vacuum limits, if we substitute k_ — > k_/£_ 
and then take the limit £_ — >• 0, we obtain the Khan-Penrose solution ^5] if <5 + = or 
the Nutku-Halil solution (201 ii S + ^ as in the previous example. The limit k_ = 0, 
5 + = leads to a known solution [S] for the collision of a step electromagnetic wave with 
an impulsive gravitational wave. The case k_ ^ 0, 5 + = corresponds to a generalization 
of the solution of Hogan, Barrabes and Bressange 



7.3 Collision of electromagnetic step-like waves of constant am- 
plitudes and polarizations 

As a final example, let us consider the collision of two step electromagnetic waves of con- 
stant amplitudes and polarizations which are not accompanied by gravitational impulses. 



OR 



It can be seen from above that, for this case, we can take 

a + (u) = cos 2 (£ + u) = cos 2 (£_t>), 

£ + (u) = -1 £_( v ) = -1 

$+(«) = e i7 + sin(£ + u) $_(t>) = e* 7 " sin(£_i>) 



(64) 



with f+(u) = 1 and /_(t>) = 1. The metric functions and additional component of the 
complex electromagnetic potential then take the forms 



H + {u) = cos 2 (£ + w), 

n+{u) = o, 



H_(v) = cos 2 (£„t>), 
$_(u) = —ie %1 ' sm{i_v) 



(65) 



where the real parameters £ ± and 7 ± determine the initial amplitude and phase respec- 
tively of the two approaching electromagnetic waves. 

In this case, the geometrically defined coordinates in the interaction region are 

a(u,v) = 1 — sin 2 (£ + w) — sin 2 (£_t>), 
(3(u, v) = sin 2 (£„t>) — sin 2 (£ + w), 

and the corresponding expressions for £ and r\ are 

£ = 1 -2sin 2 (£ + w), r] = -1 + 2 sin 2 (£_v). (66) 

Calculation of the "in-states" *Sf + (u, w) and *&_(v,w) As in the previous exam- 
ples, with the characteristic initial data ([Hljl . we can calculate, first the matrix coeffi- 
cients U(£, 77 = —1) and V(£ = 1,7/), and then the corresponding fundamental solutions 
*H? + (u,w) and *&_(v,w) of the linear systems f!26[) on the characteristics 7] = —1 (i.e. 
f = 0) and £ = 1 (i.e. w = 0). The result of these calculations for & + (u, w) and its inverse 

is _ 

*+(«, w) = i|> 0+ (u, iu) <g> k + (w) + M 0+ (u, w), k + (w) = (0, 0, 1) , 



/ ie %1+ sin(£ + u)\ 



2(w 



3*7+ 



V 



2(w - 1) 
1 



M 0+ (u,w) 




J 



where X + (u, w) is defined in (|5Uj). The inverse matrices possess the expressions 



* , 1 («,«;) = l+(w) ®<p {)+ (u,w) + N 0+ (u,w) 



'0 X 

l + (w) = -4(w 2 - 1) I 



cp 0+ (W, W) 



A+(tt,w) 
'4(w 2 - 1) 



(0, 0, 1) 



/ ze* 7+ sin(£ + u) \ 

2(w-r 

No+(u,iu) = n 1 e i7 +sin(£ +M ) 



1 
\0 



2(^-1; 





/ 



The equivalent expressions for *&_(v,w) and its inverse have the same structure except 
that 



i\> _(v,w) = A_ 1 (w,w) 



/ ie tJ - sin(£_v)\ 
2(w + l) 



V 



e n ~ sin(£^v) 

2(w + l) 
1 



N n Jv,w) 



J 








ie n ~ sm(£_v) 




2(w + l) 





1 




2(w + l) 


V> 









Construction of the integral evolution equations. We now need certain fragments 
of the structures of the matrix functions & + (u,w) and &_(v,w) on the cuts L±. These 
are the conserved monodromy data vectors 

k + H = (o, o,i), k_H = (o,o,i) 

and the jumps 4* 0+ and c)j _ of the column-vector functions \p 0+ at r + G L + and ip _ at 
r_ G L_ respectively: 



/ ie ll+ sm(£ + u)\ 



4>o+(u,t + ) = [A, ] r 



V 



2(r + - 1) 
e n+ sin(£ + u) 

2(r + - 1) 
1 



^ ie n_ sin(£_f )\ 

2(r + + l) 
e J7_ sin(£_t>) 



/ 



V 



2(r + + 1) 
1 



/ 



which are the characteristic initial values for the unknowns in our integral evolution 
equations. After that, we calculate the dynamical monodromy data vectors 

m + (v,( + ) = A_(C+)(0, 0, l), m_(u,C_) = A + (C_)(o, 0, 1 

Then, for the kernels S± defined in (|HHJ), we obtain the expressions 



S+(u,t + ,(-) 



A + KC) 

Z7r(C- -T + ) : 



m(C+ - r_) 



With these expressions, the structure of the integral evolution equations 
following dependence of the unknown functions on the spectral parameter: 



implies the 



<i> + (r + ) = [\- 1 ] T+ \Z 1 (r + ) 



'0 N 



A + 




B + 


t+ - 




r + + l 


A_ 




B 


T_ — 







*-(^) = [A: 1 kA; 1 (r_)<| |o| 



where the vector coefficients A± and B± are independent of r ± but can be functions of u 
and v. Substituting these expressions back into the linear integral equations, we obtain 



A+ = A_ = ~2 en+ sin(£ + u) coa(£_v) | i I B = B 



— -e %1 ~ cos(£ + w) sin(£_t>) ( —i 





an 



This solution of the integral evolution equations, in accordance with (|53J) leads to the 
following expressions for the Ernst potentials: 

£ = — 1, $ = e n+ sin(£ + u) cos(£_v) + e %1 ~ cos(£ + u) sm(£_v). 

In the particular case in which the electromagnetic waves possess aligned polarizations 
(i.e. for 7+ = 7_), this solution reduces to the Bell-Szekeres solution while, for 7+ 7^ 7_, 
it coincides with another solution of Griffiths (see [3]) with the metric functions 

1 - sin 2 (£ + w) - sin 2 (£_v), 

cos 2 (£ + u) cos 2 (£_v) + cos 2 (£ + u) cos 2 (£_v) — \ cos(7 + — 7_) sin(2£ + w) sin(2£_t>), 

sin(7 + — 7_) sin(2£ + w) sin(2£_t>) 
2H(u,v) ' 

8 Concluding remarks 

We have presented in this paper a general scheme for the construction of solutions of 
the characteristic initial value problem for the Einstein or Einstein-Maxwell field equa- 
tions for vacuum or electrovacuum space-times respectively with two-dimensional spatial 
symmetries. Specifically, we have demonstrated the application of this technique in the 
construction of solutions of the characteristic initial value problem for the collision and 
subsequent interaction of given plane gravitational or gravitational and electromagnetic 
waves with distinct wavefronts, which are propagating initially towards each other in a 
Minkowski background. 

The main feature of this approach is that it enables a solution to be constructed 
starting directly from given characteristic initial data. In the case of plane wave collisions, 
these data can be determined via the matching conditions from certain functions and 
parameters that are specified on the two intersecting null characteristics which correspond 
to the wavefronts of the approaching waves. Then, to construct the corresponding solution 
in the wave interaction region, it is necessary to work through three essential steps. 

• First, we must solve two systems of linear ordinary differential equations with a 
spectral parameter which determine the characteristic initial values & + (u, w) and 
ty-(v,w) of the matrix function ^f(u,v,w). 

• Next, the solutions *& + (u,w) and *&-(v,w) should be used in the construction of 
the kernels and right hand sides of the integral evolution equations derived above, 
and then we have to solve these integral equations. 

• As the last step of our construction, we have to calculate the quadratures which 
express all components of the solution of the characteristic initial value problem in 
terms of the solution of the integral evolution equations. 

It is clear, that these three steps solve our (effectively two-dimensional) characteristic 
initial value problem, at least in principle. However, it is also clear that a practical 
realization of all these steps, for more or less nontrivial initial data, may meet a number 
of technical difficulties. First of all, a problem which frequently arises is that, for given 
initial data, the ordinary linear differential equations for ^f + (u, w) and (or) *&-(v, w) may 
not possess explicit solutions. The algorithm obviously fails in this case. However, even 



a(u, v) = 
H(u,v) = 

Q(u, v) = 
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if we are successful in obtaining an explicit solution of these equations, we may meet a 
similar difficulty in the solution of the corresponding integral evolution equations. 

Despite the above-mentioned difficulties, however, the algorithm presented here can 
work in practice, as is demonstrated in the examples given above. Of course, as may be 
noted with a grain of salt, most of our examples correspond to physically important but 
well known solutions. Thus, it has not been demonstrated in the paper that the suggested 
algorithm can work effectively as a new explicit solution generating (or more correctly, - 
solution constructing) technique, although the general solution given in section 7.2 is new. 
It may be noted however, that our examples were selected using the simplest physically 
significant structures of the approaching waves. Some further physically reasonable but 
more complicated examples lead to differential and integral equations which seem not to 
admit explicit solutions in terms of elementary functions. In these cases further work may 
be useful. On the other hand, it is not difficult now to find many formally new examples 
for specially constructed, but more artificial, explicitly "integrable" initial data. However, 
looking for such examples is not of interest for us because it returns to an approach which 
is opposite to our strategy of construction of solutions starting from the initial data. 

It may thus occur that our approach can be useful for a direct construction of some 
new physically interesting (exact or possibly approximate) examples. However, it must 
be pointed out that the main purpose of our construction is a systematic description of 
the steps which are necessary to actually construct the solution from the initial data. 
In addition, this approach to the solution of the characteristic initial value problems for 
vacuum and electrovacuum fields in space-times with two commuting isometries suggests 
some general frameworks for a consideration of various questions concerning the dynam- 
ics of these fields. Among these, we can mention for example, the relation between the 
structure of the initial waves and the asymptotic behaviour of the solution in the wave 
interaction region near the singularity, the conditions for the approaching waves leading to 
the creation of an (unstable) horizon instead of a singularity etc. It may also be mentioned 
that the basic constructions of the approach suggested above, such as the spectral prob- 
lems which give rise to the characterisation of every solution in terms of its monodromy 
data and the linear integral evolution equations, can be used for a similar consideration 
(and possibly for a solution) of the analogous problems for the collision and interaction 
of non-plane waves or various waves propagating through non-Minkowskian backgrounds. 
All these questions seem to be very interesting topics for further investigations. 
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